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Abstract. Optomechanical systems are a promising candidate for the implementation of quantum 
interfaces for storing and redistributing quantum information. Here we focus on the case of a high- 
finesse optical cavity with a thin vibrating semitransparent membrane in the middle. We show that 
robust and stationary optomechanical entanglement could be achieved in the system, even in the 
presence of nonnegligible optical absorption in the membrane. We also present some prehminary 
experimental data showing radiation-pressure induced optical bistability. 
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The realization of efficient quantum networks requires the presence of quantum inter- 
faces able to store, retrieve and redistribute quantum information. Some demonstrations 
have been realized with atomic ensembles, able to store both binary [1] and continuous 
variable (CV) information [2] . Another candidate, which could be well integrated with 
both atomic and solid-state devices, is represented by optomechanical systems, in which 
information can be encoded in a vibrational degree of freedom of a resonator. 

Various cavity optomechanical systems have been recently developed [3, 4], but we 
have focused here, both theoretically and experimentally, on the setup formed by a stan- 
dard Fabry-Perot cavity with a thin semitransparent membrane placed at its center rep- 
resenting the mechanical element [5]. A given cavity mode with annihilation operator a 
([a,a^] = 1) is driven by an intense laser with frequency ft)/ and excites with its radia- 
tion pressure several vibrational modes of the membrane. However, a single mechanical 
mode can be considered when a bandpass filter in the detection scheme, centered around 
an isolated mechanical resonance, is used, and coupling between the different vibrational 
modes can be neglected. The system can be then described by two harmonic oscillators 
with Hamiltonian 

H =^{p^ + q^)+haM^)Wa + i^E{a'^e-''"" -ae'""'), (1) 

where ft),„ is the frequency of the selected mechanical mode of the membrane, q and p 
([q,p] = i) are the dimensionless position and momentum operators associated with it, 
E is related to the input laser power ^ and the mode bandwidth without the membrane 
Kq by I^I = y^l^Ko/hcoi, and (Oc[z{q)] is the frequency of the driven cavity mode. 
The latter depends upon the membrane position along the cavity axis z{q) which in 
turn depends upon the coordinate q of the membrane mode because one can write 



z{q) = Zo + -^o<?5 where zq is the membrane center-of-mass position along the cavity 
axis and xq — ^JlijnuE^, with m the effective mass of the mechanical mode [5,6]. This 
position dependence implicitly contains the radiation pressure interaction and 

can be obtained by solving the wave equation for the optical field within the cavity in 
the presence of the membrane with thickness and index of refraction hm- In the case 
of a symmetric cavity with equal mirrors and length L, and if the membrane is placed 
close to the cavity waist at the cavity center (z — 0), the frequency can be written as 



ft)^, + — arcsin 

1^ 



(-l)^cos(2^oz)(4-l) 
4nl^cot^{nMkoLa) + {nlj + 1)^ 



(2) 



where (Ot, does not depend upon z and its explicit expression is not relevant here, 
ko = COq/c is the mode wave-vector without the membrane, and p is the longitudinal 
mode index, i.e., the number of half- wavelengths within the cavity. 



Langevin equation description 

The dynamics are determined not only by the Hamiltonian of Eq. (1), but also by the 
fluctuation-dissipation processes affecting both the optical and the mechanical mode. 
The mechanical mode is affected by a viscous force with damping rate 7,,, and by a 
Brownian stochastic force with zero mean value ^ (t), whose correlation function, in the 
limit of high mechanical quality factor Qm — C),n/7m^ can be written as [7, 8] 



(3) 



where no = {exp{h(Om/kBT()} — l)~Ms the mean thermal excitation number at the mem- 
brane reservoir temperature Tq, and 5'{t — t') denotes the derivative of the Dirac delta. 

Then, due to the nonzero transmission of the cavity mirrors, the cavity field decays at 
rate Kq and is affected by the vacuum radiation input noise aQ{t), whose only nonzero 

correlation functions is given by [8] {d^{t)dQ'\t')) = d{t-t'). The presence of the 
membrane and of its non-negligible optical absorption, described by the imaginary 
part of the refraction index, hm = + inj, provides an additional loss channel for the 
cavity photons, together with the associated vacuum optical input noise af{t). This 
is a nonlinear dissipative process which affects also the mechanical mode, because 
the photon absorption rate by the membrane depends upon the resonator position q 
according to K"i[z((5r)] = Im{C()c[z(<7)]}- In particular, this process is responsible for 
an additional stochastic force on the resonator, describing membrane heating due to 
absorption. Adding all these damping and noise terms to the Heisenberg equations 
of motion associated with the Hamiltonian of Eq. (1), one gets the following set of 
nonlinear quantum Langevin equations (QLE), written in a frame rotating at ft)/. 



q = COmP, 



(4) 



(5) 



d = -/ [G)c[z{q)\ -GH]a-[KQ + /ci [z(^)]] a + E + ^/iK^d^ + V2k:i , (6) 

where denotes the derivative with respect to q. In order to reach the regime in which 
quantum effects of the optomechanical system are visible one needs a strong radiation 
pressure interaction and this requires an intense intracavity field, which is achievable 
with a large cavity finesse ^ and with enough driving power. In this case the operation 
point is given by a classical steady state characterized by a coherent intracavity field 
with amplitude tty (|«s| ^ 1) and a deformed membrane with a new stationary position 
qs, satisfying the conditions 



_ dq(Oc[z{qs)]\as? 



\0Cs\^ = ^ 7- (8) 

[Ko + Ki[z{qs)]Y + [<»/ - oMqs)]] 

Eq. (8) is a nonlinear equation determining \as\, due to the dependence of qs upon 
\as\ itself given by Eq. (7), and which may show optical bistability, i.e., the presence 
of two simultaneous stable solutions above a given threshold for the input power 
as demonstrated in [9] and also with Bose-Einstein condensates [10]. We have veri- 
fied experimentally the presence of radiation pressure-induced optical bistability in the 
membrane-in-the-middle scheme by considering a symmetric Fabry-Perot cavity with 
length L = 9 cm, waist wq = 130 jim, driven at A = 1064 nm by a Nd:YAG laser with 
input power ^ = 30 mW, and inserting at its center a macroscopic SiN membrane with 
thickness = 500 nm and side length 3 mm. Optical bistability due to radiation pres- 
sure has been observed by scanning the driving laser frequency (0/ across the resonance 
while recording the cavity transmission on a InGaAs photodiode, either from the red 
and from the blue side of the cavity resonance (see Fig. 1). An evident hysteretic cycle 
is observed, and the data are consistent with a cavity finesse ^ = 8000 and an effective 
optomechanical coupling dq(Oc[z{qs)] ~ 6 Hz. 



Quantum dynamics of the fluctuations 

The quantum behavior of the optomechanical systems could be detected in the dy- 
namics of the fluctuations around the classical steady-state described above. Rewriting 
each Heisenberg operator of Eqs. (4)-(6) as the classical steady state value plus an addi- 
tional fluctuation operator with zero mean value, and neglecting nonlinear terms in the 
equations (which is justified whenever loc^l ^ 1), one gets the following linearized QLE 
for the fluctuations 

dq = (Om^p, (9) 

5p = -[cOrn + d^COMqs)]\asn5q-rrndp + G5X + ^ + ^^^^^k^^ 
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FIGURE 1. Plot of cavity transmission for two laser frequency scans, either from the red (red curve) 
and from the blue (blue curve) side of the cavity resonance relative to a TEMqo mode. 



dX = -[KQ + Ki[ziq,)]]5X+AdY -V2dgKi[ziq,)]as5q 

+^oxt + V^T^M^s)]xr, (11) 

5Y = -[Ko + Ki[ziqs)]]dY-AdX + G5q+,/2FoYt+V^^^MQJ]Yi"- (12) 

We have chosen the phase reference of the cavity field so that a, is real and positive, 
and we have defined the detuning A = (Oc[z(<?.s)] — W/, the cavity field quadratures 
dX = {5a + 5a^)/\/2 and 5Y = {5a — 5a^)/i\/2 and the corresponding Hermitian 
input noise operators Xj" = {af + af^) 1^2 and Fj" = {a'f - a'f) /iV2, j = 0, 1 . The 
linearized QLE show that the mechanical mode is coupled to the cavity mode quadrature 
fluctuations by the effective optomechanical coupling 



G = -dqG}c[z{qs)]asV2 = -2dg(0c[z{qs)] 

which can be made very large by increasing the intracavity amplitude a^. In addition to 
cavity decay and noise terms related to K'i[z(^i)], optical absorption by the membrane 
is responsible also for two additional terms in the linearized QLE of Eqs. (9)-(12) 
which are usually neglected in standard treatments: the noise term proportional to 
Fj" describing heating of the mechanical resonator, and the dissipative coupling term 
V2dqKi[z{qs)]0Cs5q. It is therefore interesting to verify if and when absorption by the 
membrane may hinder achieving the quantum regime for the optomechanical system 
under investigation. In fact, it has been already shown that such a system may achieve 
a stationary state showing simultaneously both ground state cooling of the mechanical 
element and robust optomechanical entanglement [11, 12]. The latter could be extremely 
useful for quantum information applications because, due to its stationary nature, it is 
robust to decoherence and available at any time. 
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h(Oi ^[Ko + Ki[z{q,)]f + A'^ 



(13) 



The stationary solution of Eqs. (9)-(12) can be obtained with standard techniques and 
it is reached provided that the system is stable. Stability conditions can be derived from 
the Routh-Hurwitz criterion [12], and are slightly different from those discussed in [1 1, 
12] due to the new terms in the linearized QLE. Stability will be discussed in more detail 
elsewhere and will be assumed to be satisfied from now on. Due to the linearization and 
to the Gaussian nature of the noise operators, the steady state is a zero-mean Gaussian 
state, completely characterized by its co variance matrix (CM). The latter is given by 
the 4 X 4 matrix V with elements Vim = («z (°°) (°°) + "m (°°) "/ (°°)) /2, where Um{°°) 
is the asymptotic value of the m-th component of the vector of quadrature fluctuations 
u{t) = {5q{t),5p{t),5X{t)^5Y{t)) . The steady state CM can be determined by solving 
the Lyapunov equation [11, 12] 



AV + VA' = -Z), 



(14) 



with A the drift matrix 
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where Kj [z{qs)] = Kq+Ki [z{qs)\ is the total cavity decay rate, and D the 4 x 4 diffusion 
matrix 
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(16) 



Equation (14) is a linear equation for the CM V which can be solved, but its solution 

is very cumbersome and will not be reported here. The CM contains all the infor- 
mations about the steady state: in particular, the mean energy of the mechanical res- 
onator is given by C/ = hcOm [{^q^) + {^P^)] /2 = ^(Om [Vn -I- V22] /2 = fiOhn {n + l /I), 
where n is the effective mean vibrational number of the resonator. Obviously, in the 
absence of radiation pressure coupling it is n = hq, where hq corresponds to the ac- 
tual temperature of the environment Tq. The optomechanical coupling with the cavity 
mode can be used to engineer an effective bath of much lower temperature T <^ Tq, 
so that the mechanical resonator is cooled. Moreover from V one can also calculate 
the optomechanical entanglement between the cavity mode and the resonator at the 
steady state. We adopt as entanglement measure the logarithmic negativity Eyy, which 
is a convenient entanglement measure, easy to compute and also additive. It is defined 



as [13] £yv = max[0,-ln2r7-], where rj = 2-^/2 L{V) - [E(y)2-4dety] 



L{V) 
ofV 



1/2- 



1/2 



and 



detVi -I- dety2 — 2detyc> with Vi,V2 and Vc being the 2 x 2 sub-block matrices 



V = 



Vi Vc 
Vl V2 



(17) 




50 100 150 200 250 300 0.0 5.7 11.4 17.1 22.8 28.5 34.2 

T (K) P(mW) 

FIGURE 2. (a) Effective mean vibrational number n (blue, full line) and logarithmic negativity (red, 
dashed line) versus temperature, at fixed input power ^ = 28.5 mW; (b) versus input power at fixed 
temperature T = 1 K. Other parameters are L ~ 0.74 mm, L^i = 50 nm, m = 9 ng, (£),„/2n ~ 10 MHz, 
Q,n = 4 X 10^ 

Ref. [12] has showed that stationary ground state cooling of the mechanical resonator 
and optomechanical entanglement could be simultaneously achieved and experimental 
results in this direction have been recently attained [14]. Choosing parameters compa- 
rable to those of Ref. [6] in the linearized QLE above, we see that simultaneous cool- 
ing and entanglement is possible also in the membrane-in-the-middle scheme, because 
membrane absorption does not represent a serious limitation. This is shown in Fig. 2 
where both n and Ef^ are plotted versus temperature (a) and versus the input power (b). 
We see that, in the absence of technical limitations, ground state cooling and entangle- 
ment could be achieved even at room temperature. 
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